Characterizing the Gravitational Wave Signature from Cosmic String Cusps 
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Cosmic strings are predicted to form kinks and cusps that travel along the string at close to the 
speed of light. These disturbances are radiated away as highly beamed gravitational waves that 
produce a burst like pulse as the cone of emission sweeps past an observer. Gravitational wave 
detectors such as the Laser Interferometer Space Antenna (LISA) and the Laser Interferometer 
Gravitational wave Observatory (LIGO) will be capable of detecting these bursts for a wide class 
of string models. Such a detection would illuminate the fields of string theory, cosmology, and 
relativity. Here we develop template based Markov Chain Monte Carlo (MCMC) techniques that 
can efficiently detect and characterize the signals from cosmic string cusps. We estimate how well 
the signal parameters can be recovered by the advanced LIGO- Virgo network and the LISA detector 
using a combination of MCMC and Fisher matrix techniques. We also consider joint detections by 
the ground and space based instruments. We show that a parallel tempered MCMC approach can 
detect and characterize the signals from cosmic string cusps, and we demonstrate the utility of this 
approach on simulated data from the third round of Mock LISA Data Challenges (MLDCs). 



I. INTRODUCTION 

Cosmic strings were first predicted to form during sym- 
metry breaking in the very early universe Strings 
with G/j, ~ 10~ 6 , where /i is the string tension and G 
is Newton's constant, were studied as possible seeds for 
structure formation and the source of anisotropies in the 
cosmic microwave background (CMB) radiation Q. Ob- 
servations of the CMB combined with other cosmological 
data sets imply that G\i < 2.3 x 10 -7 , which rules out cos- 
mic strings as the source of structure in the universe [|| . 
This bound makes field theoretic cosmic strings poor can- 
didates for detection by current gravitational wave ob- 
servatories, though the next generation of detectors will 
be sensitive to gravitational waves from field theoretic 
strings Q. 

Interest in detecting gravitational wave signatures 
from cosmic string networks has increased following 
the realization that cosmic supcrstrings can be pro- 
duced in a variety of string inspired inflationary mod- 
els H i, 0, 1, & M, 03, El ■ Networks of cosmic super- 
strings can evolve quite diffently than their field theoretic 
cousins, making them more promising candidates sources 
of gravitational waves. Cosmic supcrstrings could poten- 
tially be observed via gravitational lensing, pulsar tim- 
ing, observations of the cosmic microwave background, or 
from the energy they radiate in the form of gravitational 
waves @. 

A string network consists of long, horizon sized strings 
and loops that tend to intersect and reconnect, forming 
discontinuities known as kinks. Cusps with high Lorentz 
boosts will also generically form on strings and loops. A 
cosmic string network will produce a gravitational wave 
background, as well as bursts of gravitational waves from 
cusps and kinks that stand out above the background. 
The strongest bursts come from cusps, where small por- 
tions of the string will be traveling near the speed of light, 
leading to the emission of a beam of gravitational radia- 
tion [l3l . [l4| . Current and future ground based gravita- 



tional wave detectors (LIGO-Virgo-GEO and Advanced 
LIGO) and the planned NASA/ESA space based detector 
(LISA) may be able to detect such signals and determine 
several properties of the string network [l5|, [n| [l7| ■ 

Here we consider how these cusp signals may be de- 
tected, and how well the signal parameters may be in- 
ferred using data from ground and space based gravita- 
tional wave observatories. The broad spectrum nature 
of the signals raises the possibility of joint detection in 
space and on the ground, and we investigate how this 
impacts parameter estimation. To our knowledge this is 
the first time that joint LIGO-Virgo-LISA observations 
have been considered for any gravitational wave source. 

The signal parameters that can be measured depend 
on combinations of the key physical parameters, so obser- 
vations of individual bursts are not enough to constrain 
the string model. Ultimately it will be quantities such as 
the event rate as a function of burst amplitude, and the 
power spectra of the un-resolved confusion background 
from more distant bursts and the decay of loops that will 
provide the strongest model constraints. In this paper 
we focus on individual bursts and defer the analysis of 
what can be learned from a global analysis for a later 
publication fPf\ . 

It has been pointed out that detecting cosmic strings 
with a microlensing quasar survey has been essentially 
ruled out by the low event rates (~ 10 10 quasars sources 
would have to be monitored for a year to expect a few 
events) and the long lensing periods (20-40 years) that 
arc predicted [l8|]. Still, loops of cosmic string with ten- 
sions in the range 10~ 10 < G/i < 10~ 6 are predicted to 
produce microlensing of stars in the local group of galax- 
ies [lj| . A detection of a cosmic string by a gravitational 
wave detector could thus be followed by an electromag- 
netic observation by looking in the direction determined 
by the gravitational wave observation for microlensing 
events. This possibility motivates accurate parameter es- 
timation from the gravitational wave detection to provide 
well determined sky location and orientation information 
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for the string. 

The paper is organized as follows: Section § |TT] de- 
scribes the gravitational wave signals from cosmic string 
cusps and our method for generating signal templates; 
Section § IIIII details the implementation of our Markov 
Chain Monte Carlo search technique; Section § IIVI pro- 
vides an illustration of the effectiveness of our search by 
applying it to the simulated data from the Mock LISA 
Data Challenge; and Section § [V] considers advanced 
ground based detectors and joint LIGO-Virgo-LISA ob- 
servations of cusp signals. 



II. GRAVITATIONAL WAVE SIGNATURE 

The gravitational wave signature from a cosmic string 
cusp is very simple. In the limit that the line of sight to 
the cusp is coincident with the axis of the emission cone, 
the waveform is linearly polarized and described by the 
power law 



h(t) = 2irA\t - U\ 



1/3 



(1) 



Here i» is the time when the observer sees the intensity 
peak. Equation ( [TJ is a good approximation for time 
seperations \t — 1*| short compared to the light crossing 
time, L/c of the feature that produces the cusp (2(j. The 
overall amplitude A is related to the distance to the cusp 
r, the string tension Gfi, and the characteristic length 
scale L as 
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(2) 



Since L and r are unknown, a measurement of A does 
not directly reveal the string tension. 

The signal is slightly more complicated for off-axis ob- 
servations. For small viewing angles a, defined as the 
angle between the line of sight and the axis of the emis- 
sion cone, the main effect is to round off the cusp wave- 
form (fT]) such that the power spectrum decays rapidly 
for frequencies greater than / max ~ 2/(a 3 L). Once again, 
the observable quantity, /max does not directly reveal the 
physical quantity of interest, L, since it also involves the 
unknown viewing angle a. 

The observed gravitational wave signal will also de- 
pend on how the waveform is projected onto the detector, 
which will depend on the sky location (9, <f>) and the po- 
larization angle ip. Thus, the detected signal depends on 
the parameters A, t*, / max , 4>, ip. These six measurable 
quantities depend on the two intrinsic source parame- 
ters (G/x, L) and the six extrinsic (observer dependent) 
quantities (r,a,t*,6,4>,ip). Clearly the observation of a 
single burst is insufficient to uniquely determine the in- 
trinsic source parameters. These can only be determined 
by considering the full gravitational wave signal from a 
string network, including the spectrum of the background 
and the rate distribution of the brighter bursts 



A. Template Generation 

Cosmic string cusp templates are easy to generate in 
the frequency domain as the time dependence of the de- 
tector response can be ignored: The effective duration of 
the burst is set by the lowest frequencies the gravitational 
wave detector can detect. For LISA this is / m , n ~ 10~ 5 
Hz, and for advanced LIGO the limit is f m in ~ 10 Hz, 
which leads to effective durations of 10 5 seconds and 0.1 
seconds respectively. In both cases the duration of the 
burst is short compared to the time scale over which the 
antenna pattern varies - a year for LISA and a day for 
LIGO. 

The frequency domain waveform for a burst of gravi- 
tational waves from a cosmic string cusp is a power law 
with an exponential decay beginning at the maximum 
frequency given by the viewing angle [20| . 



Hf) = 



Af-i 

Af -l e i-f/U 



f < fn 
f > fn 



(3) 



When convolved with the instrument response and 
considered relative to the instrumental noise spectrum 
the effective duration of the bursts is even shorter than 
the simple estimates based on the low frequency limit. 
Figure [1] shows that for LISA observations 99.9 + % of 
the full signal to noise ratio (SNR) is reached for obser- 
vation times of less than 400 seconds, so the stationary 
detector approximation is very accurate. 

The SNR for a template h is calculated: 



SNR = (h\h) 



1/2 



(4) 



where the noise weighted inner product for the indepen- 
dent data channels (3 over some observation time T f> s is 
defined as 
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and S n (f) is the one-sided noise spectral density in each 
channel. 



1. LISA Instrument Response 

We adopt the standard MLDC ecliptic coordinate sys- 
tem with origin at the barycenter. The individual data 
streams from the six LISA phase meters can be combined 
to cancel out the laser phase noise and form Time De- 
lay Interferometry (TDI) variables [ll[. The data sets 
for MLDC Challenge 3.4 included both the phase me- 
ter outputs and the complete set of Michelson style TDI 
variable {X, Y, Z}. The latter can be used to construct 
three noise orthogonal data streams that are similar to 
the {A, E, T} variables described in Ref. [1^. Each data 
stream s contains the response to the gravitational wave 
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FIG. 1: The SNR as a function of observation time for three 
different sets of parameter values for cosmic string cusp bursts 
of gravitational waves from the MLDC. 



FIG. 2: The low frequency approximation to the LISA re- 
sponse compared to the full response for the same source. 



signal and additive, stationary, Gaussian distributed in- 
strument noise: sg = hp + np. 

The gravitational wave response is computed by con- 
volving h(x, /) with the static limit of the LISA instru- 
ment response [23| . The various time delays in the re- 
sponse appear as phase shifts that give rise to transfer 
functions of the form 

T l3 (/) = ^inc[^-(l - k ■ nj)] exp{- l [^-(3 + k ■ 

+ ^inc[^-(l + k ■ hj)] exp{- l [^-(l + k ■ hi)]} 

where k is the propagation direction of the gravitational 
wave, fij is the unit vector pointing from spacecraft i to 
spacecraft j, and /» = l/(27rL) is the LISA transfer fre- 
quency for arm length L. The transfer function has zeros 
at frequencies that depend on the propagation direction. 
This gives LISA angular resolution for burst sources since 
bursts from different sky locations (6, <fi) will project dif- 
ferently onto the detector. 

It is the high frequency component of the gravitational 
wave signal from a cosmic string cusp that determines the 
resolution of the source on the sky. Cusp seen with large 
viewing angles result in signals with frequency cut-offs 
below the transfer frequency, and as a result, little or no 
sky resolution. Figure [5] shows the LISA instrument re- 
sponse in the low frequency limit compared to the full 
LISA response to the same source. The specific location 
of the minima of the full response depend on the sky lo- 
cation of the source and angular resolution is achieved 
because templates with different sky locations have dif- 
ferent response shapes. 



2. LIGO-Virgo Instrument Response 

The LIGO and Virgo instrument response is similar to 
the low frequency limit of the LISA response, save for a 
factor of \/3 because of the equilateral triangle configu- 
ration of LISA as opposed to the 90° orientation of the 
ground based detector arms. 

The orientation and location of each detector is most 
easily defined in an Earth fixed coordinate system with 
origin at the geocenter [24| . The Earth can be considered 
stationary for the short duration that the cusp signal is 
in the LIGO-Virgo band. This allows us to generate the 
waveform templates in the frequency domain with the 
different times of arrival at each detector appearing as 
relative phase shifts in the signal. 

During the epoch of joint LISA-LIGO- Virgo observa- 
tions the ground based detectors are expected to be op- 
crating in advanced configurations, so the baseline ad- 
vanced LIGO wideband noise curve was used for the ter- 
restrial network. 



III. MCMC TECHNIQUES 
A. Markov Chain Monte Carlo 

The Markov Chain Monte Carlo search method pro- 
vides a powerful tool for searching large parameter spaces 
without creating a large grid of templates [25|, [2(| . 

Instead of stepping along a grid to find the best fit to 
the data, a template h{x) is generated at some place in 
parameter space, x, and compared to the data to form the 
likelihood C(x). Each of the parameters is then varied by 
some amount to reach a new place in parameter space, 
y. The quality of the fit at each location is evaluated and 
the Hastings Ratio H is then formed to compare the two 
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templates (27j . 

n(f)£(f)q(y|f) W 

£(A) = e _ 3( s - /l (^)l s -' i (^)) . (7) 

The Hastings ratio is built with n(x), the prior dis- 
tribution for the template parameters, and q(x\y), the 
proposal distribution, which is the function that gener- 
ates proposals for moves from x to y. Once a jump is 
proposed, the parameters y are adopted with probabil- 
ity a = min[l,7J]. When y is a better fit to the data 
than x, H is greater than 1 and the parameters y are 
adopted. If y is a worse fit to the data, H is less than 1, 
and there is a probability equal to H that the parameters 
y will be adopted. This process is repeated until some 
convergence criteria is reached, and the parameter values 
at each iteration are stored to produce a useful map of 
the parameter space. 

A properly constructed MCMC search will sample the 
parameter space with the number of iterations spent at 
each parameter value proportional to how well that value 
fits the data. The chain of values for each parameter is 
used to form histograms that represent the marginalized 
posterior distribution functions (PDFs) that show the 
resolution and expected error for the parameter. 

The signal was parameterized by x — > 
{lnyl, t*, In / max , 0, 4>, ip}, and the priors were taken 
to be uniform in these quantities, save for 0, where a 
uniform sky distributions is given by 11(0) = |sin(0). 



The elements of the Fisher matrix are given by 

I\(A) = (MM (8) 

where hi = d\ t h. 

A Gaussian Fisher proposal distribution is good for 
local exploration of a parameter space, but more global 
proposal distributions are needed to ensure full coverage 
of the parameter space. 

2. Detector Symmetry Based Proposals 

Understanding the symmetries associated with the de- 
tection of bursts provides another effective proposal dis- 
tribution. Any burst short enough to treat the LISA 
detector as stationary will have a degeneracy such that 
sky locations related by a reflection in the plane of the 
detector with produce an identical response. There are 
additional symmetries in the low frequency limit that 
result from 120° rotations in the plane of the detector. 
These symmetries are broken at higher frequencies by the 
slightly different arrival times of the gravitational waves 
across the LISA array. Rotations thus produce two sets 
of secondary maxima (Figure [3]) . 

{d4,t, t ,cp) 

n/ 



B. Proposal Distributions 

An MCMC search is guaranteed to eventually con- 
verge on the posterior distribution regardless of choice 
of proposal distribution. In practice the performance of 
the algorithm is quite sensitive to the choice of proposal 
distribution, and care must be taken to ensure that the 
chains do not get stuck on local maxima of the PDF. We 
employed several techniques to ensure rapid exploration 
of the full parameter space: local coordinate transfor- 
mations to uncouple the parameters; moves that exploit 
symmetries of the likelihood surface to encourage jumps 
between local maxima; and parallel tempering to encour- 
age wide exploration of the posterior. 

1. Gaussian Proposals 

The Fisher Information Matrix Y can be used to pro- 
pose a jump to parameter values y based on the Gaussian 
approximation to the Likelihood at x. These proposals 
efficiently explore local maxima by calculating the eigen- 
vectors and eigenvalues from the Fisher Matrix to pro- 
pose appropriate sized jumps along independent eigen- 
directions rather than the correlated coordinate direc- 
tions. 



n \ 

(0',0X<P') 

FIG. 3: Reflection in the plane of the detector gives a degen- 
erate set of parameter values. Rotations in the plane of the 
detector give two sets of secondary maxima. 

The MCMC searches for cosmic string cusps include 
proposal jumps from x to the set of parameters y that 
give the degenerate detector response. The mapping be- 
tween these sets of parameters involves the sky location, 
polarization angle, and the time of arrival at the solar 
barycenter. The symmetry reflection and rotations are 
included as possible jumps, but the rotations are rarely 
accepted unless the chain happens to be at a secondary 
maxima and the rotation takes the chain to one of the 
two primary maxima. 

Our Fisher matrix prediction uses the LISA symme- 
try and calculates the approximation to the response at 
the two degenerate sky locations. While the likelihood is 
identical at the two locations in parameter space, the cur- 
vature of the likelihood surface is different. The Fisher 
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FIG. 4: A sky location histogram of a Parallel Tem- 
pered MCMC search for MLDC source 3.4.2 on the full 
sky with quartile contours and a Fisher matrix approx- 
imation for the source location below for comparison 
(marginalized over the other source parameters). The 
all-sky figures use the HEALPix pixelization of the sky 
( |http: / /healpix .jpl. nasa.gov ) . 




FIG. 5: A sky location histogram of a Parallel Tempered 
MCMC search for MLDC source 3.4.4 on the full sky with 
quartile contours and a Fisher matrix approximation for the 
source location below for comparison (marginalized over the 
other source parameters). 



matrix is useful for driving jumps in the MCMC searches, 
but it is not a perfect prediction of the PDF for the 
source. The maps shown in Figures [4] and [5] compare 
the marginalized PDFs for the sky location derived from 
the Fisher matrix approximation and from Markov Chain 
Monte Carlo explorations of two bursts in the Mock LISA 
Data Challenge 3.4 training data. 

We compared MCMC runs with and without the sym- 
metry based jumps, and while parallel tempering did al- 
low the chains to transition between the various maxima, 
the mixing was greatly enhanced by including the sym- 
metry based jumps. A chain with such symmetry jumps 
moves freely between the two degenerate locations in pa- 
rameter space, exploring both peaks and producing bi- 
modal histograms for each source parameter. 

The use of a variety of proposal distributions is essen- 
tial for an efficient search of the parameter space. In 
the case of cosmic string cusp gravitational wave sources 
the symmetry considerations are especially important to 
avoid getting stuck on a single maxima. The detection of 
a cosmic string cusp by LISA must include the two de- 
generate answers, with the possibility of the degeneracy 
being broken by a simultaneous detection by the ground 
based gravitational wave detector network. 



C. Parallel Tempering 

Global exploration of the parameter space is enhanced 
by creating a set of parallel chains with likelihood sur- 
faces at different "Temperatures" T such that 



(9) 



Chains that explore surfaces with T ^> 1 tend to take 
bigger steps since the contrast between maxima and min- 
ima is decreased, and this encourages wider exploration 
of the parameter space 0, [2^] . The chains can exchange 
parameters according to the Hastings ratio 



Hpt = 



£ a (x b )£ b (x a ) 

Ca{x a )Cb{xb) 



(10) 



for chains with temperature T a and Tb and parameters 
x a and Xb, respectively. 

We implement the parallel tempering method for Nq 
chains with the T values given by 



Ti = (AT)* 



where 



AT = (Tmax) ' 



(11) 



(12) 



Each chain explores its own likelihood surface until an 
exchange of parameters with a neighboring chain is pro- 
posed as a step in the MCMC algorithm. If the proposal 
is accepted, the two chains trade parameter values and 
continue to explore from this new location. Ideally, the 
chain with the highest T value has an effective likelihood 
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surface that is smooth enough that the chain makes bold 
jumps and freely explores the entire range of parameters. 
The discovery of a favorable location in parameter space 
is then propagated down to the T = 1 chain as exchange 
proposal jumps are accepted by adjacent chains. Only 
the T = 1 chain samples the true PDF and is used to 
produce the parameter histograms. 

There are several trade offs that go into the choice of 
AT and Nc- The hottest chain should have an effec- 

— 1/2 

tive signal-to-noise ratio SNR c ff = SNRT ma i < 5 to 
ensure complete coverage of the parameter space, and 
the temperature increment between chains should not be 
much larger than AT < 2 to ensure good mixing between 
chains. For example, to explore a signal with SNR = 30, 
a good choice of parameters would be AT = 1.4 and 
Nc = 12. While the cost per iteration of a parallel tem- 
pered MCMC algorithm is Nc times larger than a stan- 
dard MCMC algorithm, the improvement in the mixing, 
especially between widely separated maxima, ultimately 
leads to a far more efficient exploration of the posterior. 



IV. THE MOCK LISA DATA CHALLENGES 

To foster the development of LISA data analysis tech- 
niques, and to demonstrate mission readiness, a series 
of Mock LISA Data Challenges has been conducted [3(| • 
The simulated data sets include anticipated sources of in- 
strument noise and a wide range of astrophysical sources. 
Challenge 3.4 contains gravitational waveforms for burst 
signals from cosmic string cusps. Data analysis tech- 
niques can be used to separate the signals from the noise 
and estimate the parameters of the sources. The MLDC 
Taskforce has created both training data sets with an 
answer key for the injected source parameters and blind 
data sets with a secret answer key. The injected param- 
eter values for the blind data sets are revealed after a 
set deadline and submissions by data analysis groups are 
evaluated for accuracy [3lj . 

Challenge 3.4 is comprised of a month long data set 
(2 samples with 1 second sampling) with cosmic string 
cusp waveforms injected with a Poisson event rate of five 
events per month. The training data set includes an an- 
swer key with the number of injected sources and their 
parameter values. The training data happens to con- 
tain five sources, but there is a good chance that the 
blind data could have between two and eight sources. 
The MLDC sources have SNR uniformly distributed in 
[10, 100] and log(/ max /Hz) uniformly distributed in [-3,1]. 
This is also the first MLDC data set with non-symmetric 
instrument noise. The cusp burst sources can be found 
using a symmetric approximation for the noise (leading 
to a small systematic bias in the recovered parameters), 
or the source parameters and the individual noise levels 
can be fitted simultaneously in the search. 

The time of arrival at the solar system barycenter (i* ) 
is highly correlated with the sky location of the source. 
A search for t* leads to poor determination of the time 



of arrival of the burst due to the inherently poor reso- 
lution of the sky location. A better choice of variable 
is the time of arrival at the guiding center of the LISA 
constellation (t a ) ■ The detector time of arrival is not as 
correlated with the sky location parameters, resulting in 
better conditioned Fisher Information Matrices to drive 
the local jumps of the Markov Chain. 



A. Training Data Results 

The training data was analyzed without reference to 
the answer key so as to mimic the steps that will be 
taken to analyze the blind challenge data. The paral- 
lel tempering technique takes care of both detection and 
characterization, so the analysis does not have to be be 
broken up into distinct stages. On the other hand, run- 
ning the search on the full ~2x 10 6 seconds of data to 
find signals with duration ~ 10 3 seconds is not very effi- 
cient, so we adopted the strategy of dividing the full data 
set into 64 segments of length 32,768 seconds. Time do- 
main filters were used to limit spectral leakage, and the 
finite response of these filters meant that signals in the 
first and last ~ 10% of each segment had to be discarded. 
To ensure full coverage, a second pass was performed us- 
ing segments offset from the first by 16, 384 seconds. 

The basic parallel tempered MCMC algorithm was 
able to both detect and characterize the cosmic string 
cusp signals in each segment. The "burn-in" time for the 
chains to lock onto the signals was shortened by several 
orders of magnitude by analytically maximizing over the 
time of arrival at the detector center and the amplitude. 
These maximizations render the chains non-Mar kovian, 
and must be turned off after the burn-in is complete, and 
the samples from the burn-in must be discarded. 

A fully-fledged Baycsian MCMC analysis of the data 
would involve conditions to decide when the burn-in 
phase was complete, the ability to search for multiple 
cusp signals simultaneously, and evidence based selec- 
tion of the putative detections. While parallel tempered 
MCMC algorithms can do all of these things, we settled 
on a less sophisticated approach that could be imple- 
mented with less effort. The first stage of the analysis was 
to search each data segment using Nc = 12 chains with 
AT = 1.55 and N = 10,000 iterations. A simple SNR 
threshold was used to decide if a source had been found 
in the data segment. Triggers with SNR = (s | /z.) 1 / 2 > 8 
were recorded for further analysis (the loudest noise trig- 
gers had SNR < 6). If a trigger was found the signal was 
regressed from the data and the search repeated (in other 
words the search is sequential rather than simultaneous). 

The initial search did not fit for instrument noise levels, 
but found all five signals in the training data and recov- 
ered the source parameters to good accuracy (Table I). 
Since the segmented data is searched twice, we expect 
to find each source twice, but one trigger happened to 
fall near the boundary between segments and was thus 
discarded. The source was found on the offset pass. 
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Source 


/max (Hz) 


SNR 


A/ (Hz) 


At (sec) 


3.4.0 


2.36e-3 


53.54 


3.98e-5 


4.64 






53.99 


7.70e-5 


5.42 


3.4.1 


1.15 


21.46 


1.11 


2.54 






21.27 


2.62 


2.94 


3.4.2 


0.46 


31.07 


0.329 


0.15 


3.4.3 


1.15e-2 


73.87 


3.82e-4 


1.16 






76.58 


9.11e-4 


0.38 


3.4.4 


2.27 


14.07 


2.25 


3.44 






14.12 


4.14 


2.45 



TABLE I: The triggers produced by a search of the MLDC 3.4 
training data. All but one of the sources was detected in both 
passes through the data (one signal happened to straddle a 
data segment boundary). The injected / max values are listed, 
along with the recovered values for SNR. The difference be- 
tween the injected and recovered parameter values for f ma x 
and £a are given in the last two columns. Note that sources 
with / m ax > /* have large A/'s. These parameter errors are 
from the first stage of the filtering - the fit improves after the 
second stage. 

The initial detection of the burst could be done at sim- 
ilar cost using the template grid approach used to search 
for cosmic string bursts in LIGO data [TB|. The main 
value of the MCMC approach is that it allows us to con- 
struct the joint posterior distribution function during the 
parameter estimation stage of the analysis. 



B. Parameter Estimation 

An MCMC search for a signal in noisy data returns the 
PDF that maps out the uncertainty in each parameter. 
This can be compared to the Gaussian approximation for 
the variance in each parameter derived from the Fisher 
Information Matrix. The five distinct triggers from the 
search phase were used as starting points for the MCMC 
runs. Data segments of length 16, 384 seconds, centered 
on the time of arrival of the burst, were used in the sec- 
ond stage of the analysis. The parallel tempering routine 
used Nc = 12 chains with a maximum temperature of 
Tmax = 125, spacing AT = 1.55 and a JV = 10 6 it- 
erations. Histograms for the six source parameters for 
MLDC training source 3.4.3 are show in Figure [51 The 
sky location parameters 9 and <fi have a bimodal distri- 
bution due to the degeneracy of two sky locations from 
the detector symmetry. The time of arrival of the sig- 
nal at the detector is very well determined, matching a 
Gaussian distribution with a variance of less than one 
second. The maximum frequency cutoff can also be well 
determined for sources such as this with / max below the 
LISA transfer frequency. 

Although not much SNR is gained in the high fre- 
quency portion of the signal, the LISA transfer function 
is essential for determining the direction to the source. 



This is illustrated quite well by training data sources 3.4.3 
with SNR = 77 and / max = 0.0115 Hz. Even this source 
with high SNR has poor angular resolution (Figure [B]). 
A new search for a source with identical parameters, and 
almost identical SNR, but with / max = 1.15 Hz reveals 
the importance of the high frequency component of the 
signal (Figure [7]). 




FIG. 6: An MCMC search for an MLDC training data source 
with high SNR but a maximum frequency cutoff below the 
LISA transfer frequency. 




FIG. 7: An MCMC search for a source with the same pa- 
rameter values as Figure H3 except for a maximum frequency 
above the LISA transfer frequency. 



V. ADVANCED GROUND BASED DETECTORS 

The wide frequency range of gravitational wave sig- 
nals from cosmic string cusps allows them to be detected 
simultaneously by space based and ground based gravita- 
tional wave detectors. The current ground based detec- 
tors have not yet reported a detection of a gravitational 
wave signal, but the next generation of detectors, such as 
Advanced LIGO, will have an enhanced ability to detect 
gravitational waves, including bursts from cosmic string 
cusps. We investigated the performance of the advanced 
terrestrial network by using the MLDC training source 
parameters, with the exception of the / max parameter, 
which we set at 500 Hz to ensure that the signal will be 
in the LIGO- Virgo band. The relative orientation of the 
ground and space based detectors was fixed by setting the 
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FIG. 8: The distribution of recovered parameter values for an 
MCMC search for MLDC training source 3.4.3 using simu- 
lated LISA data. The Fisher Matrix Gaussian approximation 
is shown in blue for comparison. 
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FIG. 9: The distribution of parameter values for an MCMC 
search of simulated data using the parameters for MLDC 
training source 3.4.0 (with / ma x boosted to 500 Hz) for the ad- 
vanced ground based detector network. The Gaussian Fisher 
Matrix approximation is included in blue for comparison. 



start of the LISA observations to the autumnal equinox 
in the year 2020. 

When just considering the ground based network we 
use right ascension a and declination d to describe the sky 
location, and the polarization basis defined in Ref. [24| . 
The time of arrival of the bursts is referenced to the Earth 
geocenter. We use the reference Advanced LIGO wide- 
band noise curve with strain spectral density 

*„(/) = io-«(*— + + 

(13) 

where x = J/215 Hz. 

The ground based detector network of the 4 km LIGO 
detectors in Hanford, WA and Livingston, LA and the 
Virgo detector in Italy, if operating at advanced LIGO 
design sensitivity should respond to the MLDC sources 
with network SNRs of 24.3, 7.9, 4.3, 35.8, and 5.7 for 
sources 3.4.0 — * 3.4.4, respectively. Without the benefit 
of joint LISA observations, sources 3.4.2 and 3.4.4 may 
not be bright enough to be detected. Simulated data 
was generated for the LIGO- Virgo network and single 
chain MCMC runs of 10 6 iterations were used to pro- 
duced the parameter histograms shown in Figures [5] and 
HTJ1 The Fisher Matrix estimates for the parameter dis- 
tributions are also shown, and are seen to agree well with 
the MCMC derived PDFs. The / max parameter is poorly 
determined, resulting in a wide spread of / max values vis- 
ited by the MCMC chains. The sky location of the burst 
is determined to an accuracy of ~ 10 square degrees. 
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FIG. 10: The distribution of parameter values for an MCMC 
search of simulated data using the parameters for MLDC 
training source 3.4.3 (with / ma x boosted to 500 Hz) for the ad- 
vanced ground based detector network. The Gaussian Fisher 
Matrix approximation is included in blue for comparison. 



A. Combining Space and Ground Based Detectors 

The broad spectrum nature of signals from cosmic 
string cusps allows for joint observations between a space 
based detector such as LISA and the network of ground 
based detectors. The LISA sensitivity range is from 0.1 
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mHz to 1 Hz while LIGO is sensitive in the range 10 Hz 
to 10 kHz. A burst from a cosmic string cusp with a cut- 
off frequency above ~ 50 Hz would produce a response 
across the LIGO-Virgo-LISA network, and the extra in- 
formation from joint detection would improve estimates 
of the source parameters. One might expect that large 
Earth-LISA baseline (up to 170 seconds delay) would re- 
sult in extraordinary angular resolution, but the rela- 
tively poor determination of the time of arrival at the 
LISA guiding center (of order 1 second) diminishes the 
effect. 

For joint observations we adopted the Barycentric 
ecliptic coordinate system used to describe LISA obser- 
vations. For the time of arrival of the burst we used the 
geocenter time of arrival (f©), as this choice gives the 
smallest correlation with the sky location uncertainty. 

Table II shows the parameter estimation improvement 
gained by combining LISA with a ground based network 
detection. The gains in sensitivity can mostly be at- 
tributed to the higher SNR. This can be understood by 
looking at the off-diagonal components of the Fisher In- 
formation Matrix, and noting that even without the LISA 
contribution the parameter correlations are already quite 
small. On the other hand, the addition of the LISA data 
reduces the sky location uncertainty to less than a square 
degree, which improves the prospects of performing a mi- 
crolensing follow-up. 



MLDC 3.4.0 


SNR 


a e (rad) 


o> (rad) 


at (sec) 


LIGO- Virgo 
LISA + LIGO- Virgo 


24.3 
76.2 


5.7e-02 
2.5e-02 


3.8e-02 
9.6e-03 


1.2e-03 
5.1e-04 




MLDC 3.4.1 


SNR 


og (rad) 


o> (rad) 


at (sec) 


LIGO- Virgo 
LISA + LIGO- Virgo 


7.9 
22.9 


1.0e-02 
4.7e-03 


6.9e-02 
4.8e-02 


4.5e-04 
2.4e-04 




MLDC 3.4.3 


SNR 


erg (rad) 


o> (rad) 


a t (sec) 


LIGO- Virgo 
LISA + LIGO- Virgo 


35.8 
84.8 


1.7e-02 
5.3e-03 


2.4e-02 
2.4e-03 


3.5e-04 
5.9e-05 



TABLE II: A joint detection by LISA and an advanced ground 
based network results in improvements in parameter estima- 
tion as shown for three sources from the MLDC 3.4 training 
data. 



VI. CONCLUSION 

A matched filter analysis using parallel tempered 
Markov chain Monte Carlo techniques can both detect 
and characterize the gravitational wave signals from cos- 
mic string cusps in simulated LISA and LIGO- Virgo 
data. Cosmic string cusp waveforms have a broad enough 
frequency range that they can be detected jointly by 
space and ground based detectors, and these joint de- 
tections significantly improve the angular resolution over 
what can be achieved by LISA alone. The addition of 
LISA observations to the ground based network leads 
to a more modest improvement that may improve the 
prospects of finding electromagnetic counterparts from 
microlensing of stars in foreground galaxies. 

The analysis presented here can be improved in several 
ways. Rather than performing the search with the SNR 
of the best fit template being used as a frequentist 
detection statistic, the Bayesian evidence for the signal 
can be computed by thermodynamic integration across 
the parallel chains [jj . The instrument noise levels 
must be treated as search parameters when computing 
the evidence, and the noise model should take into 
account the non-symmetric noise levels found in the 
Challenge 3.4 data. The signal model can be extended 
to account for multiple overlapping bursts, which would 
be an improvement on the sequential regression used in 
the current analysis. A more realistic analysis should 
also take into account the confusion background that 
would accompany the bright burst signals. 
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